Parity- Violating Interaction Effects I: the Longitudinal 
Asymmetry in pp Elastic Scattering 



J. Carlson 

Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 



R. Schiavilla 

Jefferson Lab, Newport News, Virginia 23606 
and 

Department of Physics, Old Dominion University, Norfolk, Virginia 23529 



V.R. Brown 

Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 

and 

Department of Physics, University of Maryland, College Park, MD 2074-2 



B.F. Gibson 

Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 

(February 5, 2008) 



Abstract 

The proton-proton parity-violating longitudinal asymmetry is calculated in 
the lab-energy range 0-350 MeV, using a number of different, latest-generation 
strong-interaction potentials-Argonne vi$, Bonn-2000, and Nijmegen-I-in 
combination with a weak-interaction potential consisting of p- and a;-meson 
exchanges-the model known as DDH. The complete scattering problem in 
the presence of parity-conserving, including Coulomb, and parity- violating 
potentials is solved in both configuration- and momentum-space. The pre- 
dicted parity-violating asymmetries are found to be only weakly dependent 
upon the input strong-interaction potential adopted in the calculation. Values 
for the p- and w-meson weak coupling constants h v ^ and h£f are determined 
by reproducing the measured asymmetries at 13.6 MeV, 45 MeV, and 221 
MeV. 

21.30. +y, 24.80.-x, 25.40. Cm 
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I. INTRODUCTION 



A new generation of experiments have recently been completed, or are presently under 
way or in their planning phase to study the effects of parity-violating (PV) interactions in 
pp elastic scattering [0], np radiative capture and deuteron electro-disintegration || at 
low energies. There is also considerable interest in determining the extent to which PV 
interactions can affect the longitudinal asymmetry measured by the SAMPLE collaboration 
in quasi-elastic scattering of polarized electrons off the deuteron [|J, and therefore influence 
the extraction from these data (and those on the proton ||) of the nucleon's strange magnetic 
and axial form factors at a four-momentum transfer squared of 0.1 (GeV/c) 2 . 

The present is the first in a series of papers dealing with the theoretical investigation 
of PV interaction effects in two-nucleon systems: it is devoted to pp elastic scattering, and 
presents a calculation of the longitudinal asymmetry induced by PV interactions in the 
lab-energy range 0-350 MeV. 

The available experimental data on the pp longitudinal asymmetry is rather limited. 
There are two measurements at 15 MeV || and 45 MeV [0], which yielded asymmetry 
values of (—1.7 ± 0.8) x 10 -7 and (—2.3 ± 0.9) x 10~ 7 , respectively, as well as more precise 
measurements at 13.6 MeV §, 45 MeV §, and 221 MeV yielding (-0.95 ±0.15) x 10~ 7 , 
(—1.50 ± 0.23) x 10~ 7 , and (+0.84 ± 0.29) x 10~ 7 , respectively, and finally a measurement 
at 800 MeV in Ref. [|10j, which produced an asymmetry value of (+2.4 ± 1.1) x 10~ 7 . 

The theoretical (and, in fact, experimental) investigation of PV effects induced by the 
weak interaction in the pp system began with the prediction by Simonius ]ll[ that the longi- 
tudinal asymmetry would have a broad maximum at energies close to 50 MeV, and that, be- 
ing dominated by the J=0 partial waves, it would be essentially independent of the scattering 
angle. A number of theoretical studies of varying sophistication followed [p^2|-[l~4]| , culminat- 
ing in the study by Driscoll and Miller ||15|| , who used a distorted- wave Born-approximation 
(DWBA) formulation of the PV scattering amplitude in terms of exact wave functions ob- 
tained from solutions of the Schrodinger equation with Coulomb and strong interactions. In 
fact, the authors of Ref. [15] investigated the sensitivity of the calculated asymmetry to a 
number of realistic strong-interaction potentials constructed by the late 1980 's. The model 
adopted for the PV weak-interaction potential, however, was that developed by Desplan- 
ques and collaborators [TTJ|, the so-called DDH model. In the pp sector, this potential is 
parametrized in terms of p- and tu-meson exchanges, in which the PV NNp and NNuj weak 
coupling constants are calculated in a quark model approach incorporating symmetry tech- 
niques like SU(6)vi/ and current algebra requirements. Factoring in the limitations inherent 
to such an approach, however, the authors of Ref. |T(| gave rather wide ranges of uncertainty 
for these weak coupling constants. 



The present work sharpens and updates that of Ref. ||15|| . It adopts the DDH model 
for the PV weak-interaction potential, but uses the latest generation of realistic, parity- 
conserving (PC) strong-interaction potentials, the Argonne v i8 JT7j, Nijmegen I fL8| , and 



CD-Bonn . Rather than employing the DWBA scheme of Ref. [JHJ to calculate the PV 
component of the pp elastic scattering amplitude, it solves the complete scattering problem 
in the presence of these PC and PV potentials (including the Coulomb potential), in either 
configuration- or momentum-space, depending on whether the Argonne Vis and Nijmegen 
I or CD-Bonn models are used. Such an approach allows us to obtain the PC and PV 
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wave functions explicitly. While this is unnecessary for the calculation reported here-the 
DWBA estimate, along the lines of Ref. Jy|, of the PV component of the pp amplitude 
should suffice-it becomes essential for the studies of np radiative capture and deuteron 
electro-disintegration planned at a later stage. 

The remainder of the present paper is organized as follows. In Sec. [Tl] the PC and 
PV potentials used in this work are briefly described, while in Sec. a self-consistent 
treatment of the pp scattering problem is provided along with a discussion, patterned after 



that of Ref. [15[], of the Coulomb contributions to the longitudinal asymmetries measured 
in scattering and transmission experiments. In Sec. [TV] the results for the asymmetry are 
presented; in particular, their sentitivity to changes in the values of the weak coupling 
constants and/or short-range cutoffs at the strong- and weak- interaction vertices is studied. 
Finally, Sec. [V] contains some concluding remarks. 



II. PARITY-CONSERVING AND PARITY- VIOLATING POTENTIALS 

The parity-conserving (PC), strong-interaction potentials used in the present work are 
the Argonne v ls (AV18) [T7j, Nijmegen I (NIJ-I) @, and CD-Bonn (BONN) @ models. 



The AV18 and NIJ-I potentials were fitted to the Nijmegen database of 1992 P0| , |2"T|] , con- 
sisting of 1787 pp data, and both produced x 2 per datum close to one. The latest version of 
the charge-dependent Bonn potential, however, has been fit to the 1999 database, consisting 
of 2932 pp data, for which it gives a x 2 per datum of 1.01 ||19[| . The substantial increase in 
the number of pp data is due to the development of novel experimental techniques-internally 
polarized gas targets and stored, cooled beams. Indeed, using this technology, IUCF has 
produced a large number of pp spin-correlation parameters of very high precision, see for 



example Ref. p2[. It is worth noting that the AV18 potential, as an example, fits the 



post-1992 and both pre- and post-1992 pp data with x 2 's of 1.74 and 1.35, respectively [|T9 



Therefore, while the quality of their fits has deteriorated somewhat in regard to the extended 
1999-database, the AV18 and NIJ-I models can still be considered "realistic". 

These realistic potentials consist of a long-range part due to one-pion exchange (OPE), 
and a short-range part either modeled by one-boson exchange (OBE), as in the BONN and 
NIJ-I models, or parameterized in terms of suitable functions of two-pion range or shorter, 
as in the AV18 model. While these potentials are (almost) phase-equivalent, they differ in 
the treatment of non- localities. AV18 is local (in LSJ channels), while BONN and NIJ-I 
have strong non- localities. In particular, BONN has a non-local OPE component. However, 
it has been known for some time |p3[] , and recently re-emphasized in Ref. ||24|| , that the local 



and non-local OPE terms are related to each other by a unitary transformation. Therefore, 
the differences between local and non-local OPE cannot be of any consequence for the 
prediction of observables, such as binding energies or electromagnetic form factors, provided, 
of course, that three-body interactions and/or two-body currents generated by the unitary 
transformation are also included p5f . This fact has been demonstrated |26| in a calculation 
of the deuteron structure function A(q) and tensor observable T 2 o(q), based on the local 
AV18 and non-local BONN models and associated (unitarily consistent) electromagnetic 
currents. The remaining small differences between the calculated A(q) and T2o(q) are due to 
the additional short-range non-localities present in the BONN model. Therefore, provided 
that consistent calculations-in the sense above-are performed, present "realistic" potentials 



3 



will lead to very similar predictions for nuclear observables, at least to the extent that these 
are influenced predominantly by the OPE component. 

As already mentioned in Sec. Q, the form of the parity-violating (PV) weak-interaction 
potential was derived in Ref. ||T6 j— the DDH model: 

9a h p a p m 



.PV 



E 

a=p,cu 



Air 



m 



|m a (l + K a )Y'(m a r)((T 1 x cr 2 ) ■ r 



+ (ai - cr 2 ) ■ [p, Y(m a r)} + } 



(2.1) 



where the relative position and momentum are defined as r = ri — r 2 and p = (pi — P2)/2, 
respectively, [...,...], denotes the anticommutator, and m and m a are the proton and 
vector-meson (p or u) masses, respectively. Note that the first term in Eq. (|2.1| ) is usually 
written in the form of a commutator, since 

i [p , Y(m a r)]_ = m a Y'(jn a r) r , (2.2) 

where Y'{x) denotes its derivative dY(x)/dx. The Yukawa function Y(x a ), suitably modified 
by the inclusion of monopole form factors, is given by 

,2 



Y(x n ) = — 




— e 



(A a /m a )x a 




1 - 




(2.3) 



where x a = m a r. Finally, the values for the strong-interaction p- and cu-meson vector and 
tensor coupling constants g a and K a , as well as for the cutoff parameters A a , are taken from 
the BONN model [|19|1 , and are listed in Table |I|. The weak- interaction coupling constants 
h pp and h pp correspond to the following combinations of DDH parameters 

K 



h 



pi 



Pa 



(2.4) 
(2.5) 



Their values, reported in Table | in the column labeled (DDH-adj), are obtained by fitting 
the available data on the longitudinal asymmetry, see Sec. |IV|. The values corresponding to 
the "best" estimates for the h Pi and suggested in Ref. [TO are also listed in Table | in the 
column (DDH-orig). Indeed, one of the goals of the present work is to study the sensitivity 
of the calculated longitudinal asymmetry to variations in both the PV coupling constants 
and cutoff parameters. In this respect, it should also be noted that, in the limit A P =A UJ and 
ignoring the small mass difference between m p and m u , were it not for the different values of 
the tensor couplings k p and k u , the p- and cu-meson terms in v PV would collapse to a single 
term of strength proportional to g p h p p p + g^h^f. 



III. FORMALISM 

In this section we discuss the pp scattering problem in the presence of a potential v given 

by 

y = v ^ + v py + v c ; (31) 

where v PC and v PV denote the parity-conserving and parity-violating components induced 
by the strong and weak interactions, respectively, and v c is the Coulomb potential. 
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A. Partial- Wave Expansions of Scattering State, T- and S'-Matrices 

The Lippmann-Schwinger equation for the pp scattering state | p, SMg)^', where p is 



the relative momentum and SM S specifies the spin state, can be written as f27 

1 



p,SM s )^=\p,SM s )£ } + 



E - Ho - v c ± ie 



p,SM s )& 



(3.2) 



where Ho is the free Hamiltonian, v = v PC + v py , and | . . .)q are the eigenstates of Ho + v c , 



E-H -v c ) \p,SM s )g> = 



(3.3) 



with wave functions given by 

<P i p % s (r) = (r\p,SM s )£ ) 

JMjL 



F L (rj;pr) r JMj Mj 

/j LSM S y£ > )\ ~yLs.j\ r ) ■ 



pr 



(3.4) 



Here Fl{ti; p) denotes the regular Coulomb wave function [28|, while the parameter r\ and 
Coulomb phase-shift <jl are given by 



77 = afi/p , 
a L = arg [T(L + 1 + ir/)] 



(3.5) 
(3.6) 



where a is the fine structure constant and /1 is the reduced mass. Finally, the following 
definitions have also been introduced: 



*3& a (p) = Y,( LM l, SM S I JM,) Y LML (p) 

M L 



(3.7) 



__ 1 r 
€LS= 2 



l + (-l 



(3.8) 



The factor 6^5 ensures that the wave functions 0^ are properly antisymmetrized. Note 
that in the limit i]=0, equivalent to ignoring the Coulomb potential, the latter reduce to 
(antisymmetrized) plane waves, 



1 
72 



3 ip-r + r\s e -ip-v 



Xmo 



(3.9) 



The T-matrix corresponding to the potential v + v c can be expressed as |27 

T(p', S"M£; p, SM S ) = T(p', S"M£; p, SM S ) + T c (p', S"M£; p, SM S ) , (3.10) 



where T c is the (known) T-matrix corresponding only to the Coulomb potential f27 |, and 
T(p', S"M£; P , SM S ) = { c \p', S'M' S \T\p, SM 3 )™ . (3.11) 



Insertion of the complete set of states | p, SMs)q into the right-hand-side of the Lippmann- 
Schwinger equation leads to 



p,sm s )m=\p,sm s )^+ y: 



dp' 1 



S'M' 



^ T(p',S'M' s ;p,SM s ) 



(3.12) 



from which the partial wave expansion of the scattering state is easily obtained by first 
noting that the potential, and hence the T-matrix, can be expanded as 



{ c\p', S'M* | v | p,SM s )^ = 2(47r) 2 ]T Y,*l>s> e LS e^' Z J f s f M ,(p') 



v(+) 



JMj LL> 



with 



J / / \ • L—L' f j Fl'{V'iP r )-v;Mjt / \i)Mr 

K'S'MP iP) = l J dr ~r r yus'jv{r)yLSj 



F L (r];pr) 



p'r 



pr 



(3.13) 



(3.14) 



After insertion of the corresponding expansion for the T-matrix into Eq. (|3.12 ) and a number 
of standard manipulations, the scattering-state wave function can be written as 



/(+) / \ a fKST^ ST^ -L' icrr \ryJMj I - \1* W L I ' S' ,Ls( r 'l P) 



(3.15) 



with 



w 



L'S 



,Ls( r ',P) 



JMj LL'S' 



8l,L'$S,S' 



FjArnpr) 
pr 



+ 



2 f°° j / ,2 F L'(v,P'r) 



TC JO 



dpp 



p'r E-p' 2 /(2n) +ie 



TZ T L'S',Lsip'\P) 



(3.16) 



The (complex) radial wave function w(r) behaves in the asymptotic region r —>■ oo as 



wi,Jr-p) 



5 a ^h { l){ri;pr) + tiff (r);pr)S£, (p) 



(3.17) 



where the label a (a') stands for the set of quantum numbers LS (L'S'), the on-shell (p' = p) 
5-matrix has been introduced, 



s I>Ap) = <W - & lipTi Jp;p) , 



(3.18) 



and the functions hS 1,2 >(ri\ p) are defined in terms of the regular and irregular (Gl) Coulomb 
functions as 



h ( l (v,p) 



(3.19) 



Again, in the limit rj = 0, F L (r);p)/p -> j L (p) and G L (r);p)/p -> -n L (p), where j L (p) 
and nii^p) are the spherical Bessel functions, and the familiar expressions for the partial 
wave-expansion of the scattering state, S- and T-matrices are recovered [^7] . 
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B. Schrodinger Equation, Phase-Shifts, and Mixing Angles 



The coupled-channel Schrodinger equations for the radial wave functions w(r) read: 

(3.20) 



dr 2 



p r 



with 



v J a ,Jr) = i L - L '2 f i jdny^v(r)y t 



Mj 
aj i 



(3.21) 



where, because of time reversal invariance, the matrix v 3 a , a can be shown to be real and 
symmetric (this is the reason for the somewhat unconventional phase factor in Eq. ( |3.21| ); in 
order to mantain symmetry for both the v PC - and v PV -matrices, and hence the S'-matrix, the 
states used here differ by a factor i L from those usually used in nucleon-nucleon scattering 
analyses). The asymptotic behavior of the w(r)'s is given in Eq. ( |3.17 ). 

The Pauli principle requires that there be a single channel when J is odd, and three 
coupled channels when J is even, with the exception of J=0 in which case there are only two 
coupled channels, 1 So and 3 Po- The situation is summarized in Table |l[ Again because of 
the invariance under time- inversion transformations of v PC + f PV , the S-matrix is symmetric 
(apart from also being unitary), and can therefore be written, for the coupled channels having 



J even, as p7 | 



s J = U T S^U , 

where U is a real orthogonal matrix, and Sp is a diagonal matrix of the form 



J D:a',a 



2iS'> 



Here 5^ is the (real) phase-shift in channel a, which is function of the energy p 
The mixing matrix U can be written as 



U 



[/(12) 

n u (ij) 

1<«<7<3 



J = , 

J > 2 with J even 



(3.22) 



(3.23) 

(3.24) 
(3.25) 



where is the 2 x 2 or 3 x 3 orthogonal matrix, that includes the coupling between 
channels i and j only, for example 





cos e{ 3 





sin e( 3 ' 










1 " 


[/(13) = 





1 





— L fc 13 













. -sin ef 3 





cos e{ 3 




-1 









Thus, for J=0 there are two phase- shifts and a mixing angle, while for J even > 2 there 
are three phases and three mixing angles. Of course, since |v PV | <C \v \, the mixing angles 
e{j induced by f PV are <C 1, a fact already exploited in the last expression above for U. 
Given the channel ordering in Table [H]> Table |T| specifies which of the channel mixings are 
induced by v PC and which by v PV . 
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The reality of the potential matrix elements v 3 , a (r) makes it possible to construct real 
solutions of the Schrodinger equation Q3.20| ). The problem is reduced to determining the 
relation between these solutions and the complex iw(r)'s functions. Using Eq. ( |3.22| ) and 
U T U = 1, the w(r)'s can be expressed in the asymptotic region as 



2_^{U ja'pe f> Up a 

P 1 

l^( u V/? e — 



U, 



p 



pr 



pa , 



(3.26) 



where the el'S' has been dropped for simplicity. The expression above is real apart from 



the exp(i<5§) 
introduced 



To eliminate this factor, the following linear combinations of the w(r)'s are 



u. 



1 T m J 

t',a _ e~ i(? /3 a '<P (jjT^ 



Pa 



{TP 



cos 5 J a Fy (tj; pr) + sin 5% Gl'(tj;pt) 
pr 



(3.27) 



and the -u(r)'s are then the sought real solutions of Eq. ( |3.20| ). 

The asymptotic behavior of the w(r)'s can now be read off from Eq. (|3.27| ) once the 
£7-matrices above have been constructed. The latter can be written as, up to linear terms 
in the "small" mixing angles induced by t> PV , 



U 



12 



J = 



12 



u 



cose{ 2 sine{ 2 e{ 3 cose{ 2 + e^ 3 sine{ 2 
-sme{ 2 cose{ 2 -e{ 3 sme{ 2 + e J 23 cose J l2 



13 



— e 



28 



J > 2, J even 



Inverting the first line of Eq. (|3.27|) , 



^^e 1 ^^ , (3.28) 
r ~ r 



and inserting the resulting expressions into Eq. ( |3.20|) leads to the (in general, coupled- 
channel) Schrodinger equations satisfied by the (real) functions uir). They are identical to 
those of Eq. ( |3.20| ), but for the u>(r)'s being replaced by the «(r)'s. These equations are then 
solved by standard numerical techniques. Note that: i) v 3 a i a = v^^, since the diagonal 
matrix elements of v PV vanish because of parity selection rules; ii) in the coupled equations 
with J even, terms of the type r v^ f « {r)ui a (r) j r involving the product of a parity- violating 
potential matrix element with a v PV -induced wave function are neglected. 
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C. Amplitudes, Cross Sections, and the Parity- Violating Asymmetry 



The amplitude for pp elastic scattering from an initial state with spin projections mi, 
m 2 to a final state with spin projections m[, m' 2 is given by 

{m x m 2 | M | vn x m^ = E (^77^, ^m' 2 | S'M' s )(^m 1 , }-m 2 \ SM S ) 

S'M' S ,SM S 11 11 

M S 'M' S ,SM S {E,9) , (3.29) 

where the amplitude M is related to the T-matrix defined in Eq. ( |3.10| ) via 

Ms'M' s ,SM s (E,e) = -^-T{p',S'M' s ;pz,SM s ) . (3.30) 

Note that the direction of the initial momentum p has been taken to define the spin quan- 
tization axis (the z-axis), 9 is the angle between p and p', the direction of the final momen- 
tum, and the energy E = p 2 /(2fi) (= p' 2 /(2fi)). The amplitude M is split into two terms, 
M = M + M c , as in Eq. ( |3.10| ). Using the expansion of the T-matrix, Eq. ( |3.13| ) with 
V L'S' ls replaced by T[, s , LS , and the relation between the S- and T-matrices, Eq. ( |3.18|) , 
the amplitude induced by v PC + v PV can be expressed as 

M s , M ^ Ms (E, 6) = V^Y. V2L + le L/S ' e LS (L'(M S - M' s ), S'M' S \ JM S ) 

JLL> 

(L0, SM S \ JM s )Y L , {Ms _ M , s) {e) S ™*s(p)- <W 6 S ,, S ^ 

6 ip 

while the partial-wave expansion of the amplitude associated with the Coulomb potential 
reads [j27f : 



p2iox _ i 

(E,e)=5 s ^sS M ' s ,Ms^J2^ r +^ e LsY L o(e) : . (3.32) 

L l P 

The differential cross section for scattering of a proton with initial polarization mi is 
then given by 

W mi {E,6) = ^E E I K m 2 I M I mim 2 ) | 2 , (3.33) 

■mi m ' im > 2 

and the longitudinal asymmetry is defined as 

-r (F B ,_ <r+{E,6)-a-{E,e) 

A{E ' 6) -a + (E,9) + a4E,9) ' (3 " 34) 

where ± denote the initial polarizations ±1/2. Carrying out the spin sums leads to the 
following expression for the asymmetry: 





'M s , M ^{E,o)W s , M , s ^{E,e) 


+ c.c. 


Y,S'M' S Y,SM S 


Ms>m' s ,sm s (E, 0) 


2 



A (E, 9) = , — s ' " „ ' , 9 i , (3.35) 



from which it is clear that the numerator would vanish in the absence of parity-violating 
interactions, since v PC + v c , in contrast to v PV , cannot change the total spin S of the pp 
pair. 

Parity-violating scattering experiments typically measure the asymmetry weighted over 
a range [9i, # 2 ] of scattering angles, 



A(E)) 



k<e<e 2 dQa(E,9)A(E, 
f 01 < e <e 2 dna(E,9) 



(3.36) 



where a = (a + + cf_) /2 is the spin-averaged differential cross section. In contrast, transmis- 
sion experiments measure the transmission of a polarized proton beam through a target. A 
cross section is then inferred from the transmission measurement. Beam particles elastically 
scattered by angles greater than some small critical angle 8q are removed from the beam, 
thus reducing the observed transmission and adding to the inferred cross section. Beam 
particles scattered at angles smaller than #0 are not distinguished from the beam and do 
not contribute to the cross section. To derive an expression for the asymmetry in this case, 
one needs to carefully consider the Coulomb contribution to the cross section-a divergent 
quantity in the limit 8q — > 0. To this end, following Ref. [25], one first defines the differential 
cross sections 



a S'M'SM s 



and 



and hence 



= I M SIM > s>SMs (E, i 
= I M S i M > stSMs (E,( 



T S'M>SM s ( E i 



12 
|2 



^S'M'SM S (E, ' 



-2 3? 



Ms>M' s ,SM s ( E i 9) M S'*M'SM s ( E i 



Ms>M'SM s (Ei 1 



+ a S'M',SM s {Ei 



&S'M' S ,SM S {.E,6) = (JS'M' S ,SM S {E 

In transmission experiments, the quantity of interest is 

/•tt/2 

vsm s ,>{E) = 2tt d6sm8 °S'M' s ,sm s { e 
Je ° S'M' S 

= CTSM S ,>( E ) + a SM s ,>( E ) > 

where &s Ms;> (E) is explicitly given by 

1 1 r-^ s 



(3.37) 
(3.38) 
(3.39) 



a 



SM S ,> 



V 
J? 



2 r 



sin 2 fl /2 cos 2 # /2 



V 



sin [277lgtg(0 o /2)] 



(3.40) 



(3.41) 



To evaluate <jsm s ,>( e ), one writes, following Ref. J29fl: 



S'M' S I 



&sm s ,>( e ) = cr SMs (E) -27t/ d6>sin6> 22 

+ 2 $t\M s , M , >SMs (E, B)M^ m , SMs {E 



M s 'M'SM S ( E , 



(3.42) 
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Application of the optical theorem to the total cross sections a and a c allows one to deduce 

47T, 

V 



a SMs (E) = a SMs (E) - a^ Ms (E) = — 3 



M : 



SM S ,SM 3 



(£7,0) 



(3.43) 



and the determination of the cross section <t S m s ,>(.E) is reduced to evaluating the integral 
on the right-hand-side of Eq. (|3.42|). For sufficiently small 8q and by appropriately taking 



the limit e —>■ in the integral of the interference term MM C *, which essentially entails 
taking the limit term by term in the partial- wave expansion of M c , one finds 



<?sm s ,>{.E) 



Air 
V 



M 



SM S ,SM S 



(E,0)e 



2ifalgsin(0o/2)-<To] 



(3.44) 



neglecting terms of order 8q and higher. Using Eqs. (|3.41 ) and (|3.44 ), the longitudinal 
asymmetry measured in transmission experiments is obtained as 



ME) 



3 


[M 10fi0 {E,0) + M OOAO {E,0)] 




3 


e^EsM s M SMs ,SM s (E,0)_ 


+ (p/^)Esm s ^ Msi> (E) 



(3.45) 



with exp(i0) = exp 2i[r] lg sin(6 l o/2) — (x ] 



D. Momentum-Space Formulation 



In order to consider the (parity-conserving) Bonn potential [[HJ , it is necessary to develop 
techniques to treat the pp scattering problem in momentum space. A method first proposed 
in Ref. |30| and most recently applied in Ref. |19[] is used here. It consists in separating the 
potential into short- and long-range parts Us and v^, respectively, 



vs + v L , 



(3.46) 



where 



v s 



6{R-r) 



v L = v c 6{r - R) , 



(3.47) 
(3.48) 



and 8(x) is the Heaviside step function, 6(x)=l if x > 0, =0 otherwise. The radius R is 
chosen large enough, so that t> PC + v PV vanishes for r > R (in the present work, R = 20 fm). 

Since Us is of finite range, standard momentum-space techniques can now be used to 
solve for the .fTg-matrix in the J channel(s): 



v J 



S;a',a 
O roo 

7T JO 



(p'; p) 

dkk 2 ^Vs; Q ', p{p'\ k ) 



V 



p 2 /(2/i) - fcV(2/i) 



Ki^ a {k;p) 



(3.49) 



where V denotes a principal-value integration, and the momentum-space matrix elements 
of the potential Us are defined as in Eq. (|3.21|) , but for the replacements v — > 11$ and 
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Fl(t}',x)/x — > jh{x)- Note that performing the Bessel transforms of a Coulomb potential 
truncated at r = R poses no numerical problem. The integral equations ( 3.49Q are dis- 
cretized, and the resulting systems of linear equations are solved by direct numerical inver- 
sion. The principal- value integration is eliminated by a standard subtraction technique [|3T 
The asymptotic wave functions associated with vs have the form: 



^S;a',a( r ;P) 



where 



h { L 2 \p) = j L {p) ±in L (p) 



(3.50) 



(3.51) 



jl and rix being the regular and irregular spherical Bessel functions, respectively, and the 
constants a a can only depend upon the entrance channel a, see the Schrodinger equa- 
tions ( |3.20D . These wave functions should match smoothly, at r=R, those associated with 
the full potential vs + Ul, which behave asymptotically as in Eq. ( [3.171) . Carrying out 
the matching for the functions and their first derivatives leads to a relation between the 
S-matrices Ss ;a ',a an d S a / tQ , corresponding to vs and vs + vl, respectively. In terms of 
i^-matrices, related to the corresponding S'-matrices via 



S J (p)= 1 + 2i fip K J (p;p) 1 - 2i up K J (p;p) 



(3.52) 



and similarly Ss and K$, this relation reads in matrix notation |32 
2fxpK J 



n -l 



G 



3 + 2fj,pN-Ki y + 2fipN'-Ki ■ G' 



x 



J + 2/xpN • Ki] \j' + 2/jpN' ■ Ki] -F'-F 



(3.53) 



where the dependence upon p is understood, and the diagonal matrices X and X' have been 
defined as 



X a ' )Q = 5 a ', a X a (R) , 
"dX Q (r) 



X' / „ = 5 n i 



dr 



(3.54) 
(3.55) 



J r=R 



with the functions 



X a (R) = j L (pR), , n L {pR), and 



pR 



pR 



(3.56) 



when X=J, F, N, and G, respectively. Once the if-matrices in the various channels have 
been determined, the corresponding S'-matrices are obtained from Eq. ( p.52|) , from which 
the amplitude Ms>m' ,sm s (E, 9) (|3.31| ) is constructed. 
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E. Matrix elements of w PV in channel J 



To evaluate the radial functions ?;^ v (r) of the PV potential in Eq. ( |3.21|) -those associ- 
ated with the PC potential are well known-one needs the matrix elements of (eri x <x 2 ) ■ r 
and (eri — <x 2 ) ■ p between spin-angle functions. Using the notation 



(J; L', S'\0\J;L,S)= J dfi y%$j O(r) y 



Mj 
LSJ J 



(3.57) 



and writing 



[<Ti - cr 2 ) ■ p = -i(<Ti - (T 2 ) 



,d Id 

r 1 

dr r <90 



(3.58) 



where the d symbol indicates that the partial derivatives must act to the right, one finds 
that the non-vanishing matrix elements are: 



(J; J, I (cr 1 x <x 2 ) • r I J; J =F 1, 1) = ±i 



(J;J,0| (0-1-0-2) - f I J; J +1,1) = ±1 



/ 2 J + 1 =f! 
/ J +1/2 

'2 J + 1 =f 1 
J+ 1/2 : 



(3.59) 
(3.60) 



J ' 1 " CT2 ) ' an 1 ' ^ ' ^ = 2 — V ^ + 1/2 



1 I (o : - o 2 ) • — I J; J, 0) = ^ J+1/2 



(3.61) 
(3.62) 



Note that the operators in Eqs. (|3.59| )- (|3.60|) are Hermitian, while those in Eqs. (|3.61|) - (|3.62|) 
are not. The complete Hamiltonian is, of course, Hermitian. 



IV. RESULTS AND DISCUSSION 

In this section we present results for the longitudinal asymmetry in the lab energy range 
0-350 MeV. The calculations use any of the modern strong-interaction potentials, either 
Argonne u 18 (AV18) g7|, or CD-Bonn (BONN) fTJ, or Nijmegen-I (NIJ-I) in combi- 



nation with the DDH weak-interaction potential parameterized in terms of p- and a;-meson 
exchanges The values for the p- and tu-meson coupling constants and cutoff parameters 



are listed in Table Q. The strong-interaction coupling constants and cutoff parameters are 
taken from the BONN potential, while the weak-interaction coupling constants ftP^ and h v J- 
have been determined by an AV18-based fit to the observed asymmetry. In Table [I] we also 
list the hP p and h v J values corresponding to the "best" estimates for the h Pi and h Wi suggested 



in Ref. [16], column labeled DDH-orig. 

The data points for the longitudinal asymmetry at 13.6 MeV, 45 MeV, and 221 MeV are 
those reported in Ref. 0,0, and their values are (-0.97±0.20) x 10" 7 , (— 1.53±0.21) x 10" 7 , 



and (+0.84 ±0.34) x 10 7 , respectively. The first point at 13.6 MeV has been obtained |33 
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by taking the weighted mean-and accounting for the square-root energy dependence-of the 
latest result from the the Bonn experiment at 13.6 MeV, as reported by Eversheim (Ref. [14] 
in Ref. pi]), and the 15 MeV result from Ref. ||. The point at 45 MeV has also been 
obtained |33| by combining results from measurements at 45 MeV || , 46 MeV, and 47 MeV 



(these last two both from Ref. [j3"3"]). The last point at 221 MeV is that reported in Ref. [|TJ. 
Finally, the errors include both statistical and systematic errors added in quadrature. 

The total longitudinal asymmetry, shown in Fig. [I] for a number of combinations of 
strong- and weak-interaction potentials, is defined as 

s[m 10i00 ( j e,o) + m 00i10 ( j e,o)1 

ME) = -tr- — -1 J , (4.1) 



J2sm s M S m s ,sm s (E, 0) 



where the amplitudes M s >m' s ,sm s (.E,9) are those given in Eq. ( |3.31|) . The expression above 
for A(E) ignores the contribution of the Coulomb amplitude, Eq. ( [3.32j ), divergent in the 
limit 9=0, and for this reason A(E) will be referred to as the "nuclear" asymmetry. Of 
course, one should note that Coulomb potential effects enter into A(E) explicitly through 
the Coulomb phase shifts, present in the partial-wave expansion for M S 'm' s ,sm s (E,9), and 
implicitly through the wave functions, from which the S'-matrix elements are calculated. The 
effect of including explicitly the amplitude induced by the Coulomb potential is discussed 
below. 

The calculated nuclear asymmetries in Fig. [I] were obtained by retaining in the partial- 
wave expansion for M s >m' s ,sm s (E, 9) all channels with J up to J max = 8. The curves labeled 
AV18, BONN, and NIJ-I all use the DDH potential with the coupling constants h pp and 
h pp determined by a rough fit to data (the AV18 is used in the fitting procedure). There 
is very little sensitivity to the input strong-interaction potential, certainly much less than 



one would infer from Fig. 1 of Ref. This is undoubtedly a consequence of the more 

extended pp and pn scattering data-base to which present potentials are fitted, as well as of 
the much higher accuracy achieved in these fits. An analysis of the extracted h pp and h pp 
coupling constants and their errors is presented later in this section. 

We also show in Fig. [I] the AV18 results which correspond to a DDH potential using the 
"best" estimates for the h pp and hff coupling constants |Tj| (values in column DDH-orig in 
Table |) , with the remaining p- and cu-meson strong-interaction coupling constants and cutoff 
parameters as given in Table Q. A number of comments are now in order. The data point at 



221 MeV essentially determines the value of h pp '. As pointed out by Simonius [34] (see also 
below), the dominant contributions to the total asymmetry in the energy range under con- 
sideration here are those associated with the 1 So- 3 Po and 3 P 2 - 1 D 2 partial waves. At energies 
close to 225 MeV the 1 So- 3 Po contribution, which can easily be shown to be proportional 
to cos ^(E; 1 S ) + (Ti(E) + a (E)} - cos [5(E; 3 P ) + a^E) + a (E)} using Eq. (|33T|) (here 
S^E; 1 S ) and 5(EfP ) are the strong-interaction phases), vanishes. As a result, the total 
asymmetry in this energy region is almost entirely due to the 3 P2- 1 D 2 contribution, which 



is known [33] to be approximately proportional to the following combination of coupling 
constants, h pp g p K p + h^g^K^. In the BONN model, the a>-meson tensor coupling constant is 
taken to be zero, and hence the data point at 221 MeV fixes h pp (for given g p , k p , and A p ). 
This is the reason for the ~ 44 % increase (in magnitude) of h p p p with respect to the DDH 
"best" estimate. 
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Below 50 MeV, however, the calculated total asymmetry is dominated by the 1 So- 3 Po 
contribution, approximately proportional to |34| h pp g p {2 + k p ) + h P fg UJ (2 + k u ). The increase 
in magnitude of h pp required to fit the point at 221 MeV, now leads to a total asymmetry 
below 50 MeV, which is too large when compared to experiment. Thus, in order 
to reproduce the 13.6 MeV and 45 MeV data points, the overall strength of the coupling 
constant combination above needs to be reduced significantly. Since g p {2 + k p ) and g UJ (2 + 
k u ) = 2g w have the same sign, this requires making the sign of h pp opposite to that of h pp . 

It is worth pointing out, though, that the changes in value for h pp and h pp advocated here 
are still compatible with the "reasonable" ranges for the h pp and h^f., determined in Ref. JIB 



Finally, we show in Fig. [L] the total nuclear asymmetry obtained in a calculation based 
on the old Reid soft-core potential [|3~5"l , [3~6l and a DDH potential using the following coupling- 
constant and cutoff values: g 2 p /^'n = 0.95, g^jkx = 20, k p = 6.1, = 0, A p = 1.3 GeV/c, 
A w = 1.5 GeV/c (these are all from the old r-space version of the Bonn potential |07D, and 
the "best" estimates for h pp and h pp . These model interactions are essentially identical to 



those employed by Driscoll and Miller in Ref. [15]. Indeed, our calculated total asymmetry 



is close to that obtained by these authors. It should be stressed that in Ref. |Tj| the strong- 
interaction phases and mixing angles were taken from Arndt's analysis of nucleon-nucleon 
scattering data |38|] rather than calculated from the Reid soft-core potential, as done here. 
This is presumably the origin of the remaining small differences between their results and 
ours. 

Figure |2| shows the total nuclear asymmetries obtained by including only the J=0 channel 
( 1 So- 3 Po) and, in addition, the J =2 channels ( 3 P2- 1 D2 and 3 F2- 1 D2), and finally all (even) J- 
channels up to J m ax=8. We re-emphasize that in the energy range 0-350 MeV the asymmetry 
is dominated by the J =0 and 2 contributions (among the latter, specifically those from the 
3 P2- 1 D 2 partial waves). 

Figure |3] illustrates the sensitivity of the total nuclear asymmetry to modifications of the 
A p and A w cutoff parameters in the DDH potential. Both cot-offs are multiplied by -R cu t, in 
each case h pp and h pp are then readjusted to approximately reproduce the AV18+DDH-adj 
combination. In the near point-like limit (i? cut =10), the asymmetry increases in magnitude 
by roughly a factor of 2 prior to adjustment. The resulting couplings used for this case are 
h pp =-14:.33 and h pp =+3.95. Results are also shown for i? cut =1.5 and 0.8, the latter is an 
extreme case where the cutoff parameters are near the meson masses used to determine the 
ranges. Nevertheless, the energy dependence of the asymmetry is in all cases very similar. 
The couplings used to generate the two other curves are: for i? cut =1.5, h pp =-15.32 and 
h pp =+3.92; and for ^=0.8, h pp =-10Q.7 and ^=+14.63. 

Figures |]-[] illustrate the effects of Coulomb contributions on the longitudinal asymmetry. 
Figure [| compares the total nuclear asymmetry defined above (curve labeled "C" ) with the 
total asymmetry obtained by ignoring the Coulomb potential altogether (curve labeled "no- 



C"). As already pointed out in Ref. [[TJ], Coulomb contributions to these (non-physical) 
quantities are rather small. 

Figure |5] compares the angular distribution of the (physical) longitudinal asymmetry 
obtained from the amplitudes M = M+M°, see Eq. Q3.35|) , with that calculated by replacing 
M — > M in Eq.( |3.35| ), namely ignoring the contribution of the Coulomb amplitude M c . The 
latter dominates at small scattering angles, and leads to the peculiar small angle behavior 
of the angular distribution shown in Fig. |5|, namely its changing of sign at small 8, and 
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its vanishing at #=0 (also observed in Ref. ||15|| ). It is interesting to note that the angular 
distribution at 230 MeV obtained by the authors of Ref. [15] is significantly different from 
that shown in Fig. | at 221 MeV. 

Figure |6| illustrates the effects of Coulomb contributions on the total longitudinal asym- 
metry measured in transmission experiments, Eq. ( |3.45| ), for various choices of the critical 
angle 6q (#o=2°, 5°, and 10°), see Eq. (|3.40|) . Coulomb contributions are substantial, par- 
ticularly for small 9q and energies below 100 MeV. We also demonstrate, in Fig. [7|, that the 
angular distributions of the longitudinal asymmetry are only weakly affected by different 
input (strong-interaction) potentials. 

Finally, we present an analysis of the extracted coupling constants and their experimental 
errors. This analysis employs the AV18 model and the BONN-derived strong interaction 
couplings and cut-offs (Table |) in the DDH potential. Given the weak sensitivity described 
earlier, only small changes should be expected with the use of other recent strong interaction 
potentials. The experimental data at low energies have been combined into the two data 
points shown in Fig. [I] at 13.6 and 45 MeV. The asymmetries are -0.97 ± 0.2 and -1.53 ± 
0.21, respectively, combining statistical and systematic errors. We also include the recent 
TRIUMF result of +0.84 ± 0.34 at 221 MeV (all in units of 10~ 7 ). 

Figure |8| shows contours of constant total y 2 at levels of 1 through 5 versus the coupling 
constants h pp and K pp . As is apparent in the figure, there is a rather narrow band of 
acceptable values for h pp and h pp at total y 2 = 1 for the 3 experimental data points. At 
this level, h pp can range from approximately -14 to -28, with a simultaneous (and strongly- 
correlated) variation in h pp from rougly -2 to +10, all in units of 10~ 7 . 



V. CONCLUSIONS 

We have performed an analysis of the pp parity-violating (PV) longitudinal asymmetry 
using combinations of modern-day strong interaction potentials and the DDH PV potential. 
The new experimental results from TRIUMF at 221 MeV, in combination with previous 
results at lower energy, provide a strong constraint on allowable linear combinations of p and 
uj PV coupling constants. Combining the statistical and systematic errors in quadrature, h pp 
is constrained by present data to approximately 35%, at the level of one standard deviation. 

The prime motivation for the present work is to initiate a systematic and consistent study 
of many PV observables in the few-nucleon sector, where accurate microscopic calculations 
are feasible. Recent measurements of nuclear anapole moments in atomic PV experiments 
are difficult to reconcile with earlier PV experiments in light nuclei using the simple DDH- 
orig model 133]. Up to now, this approach has been extremely limited by the available 



data. 

The present experimental data in the few-nucleon sector remain rather sparse, with 
primary constraints coming from the pp longitudinal asymmetry analyzed here and mea- 
surements of the PV asymmetry in p — a elastic scattering ||40|| . A variety of new results are 



expected in the next few years, however, including 1 H(n, j) 2 !! 0, the neutron spin rotation 
in Helium |j41 |, and electron scattering measurements at Bates and JLAB. The combination 



of these diverse experiments should finally yield a coherent picture of the NN PV interaction 
at the hadronic scale. 



16 



ACKNOWLEDGMENTS 



We wish to thank R. Machleidt for providing us with computer codes generating the 
latest version of the Bonn potential, and for illuminating correspondence in regard to the 
treatment of the Coulomb potential in momentum space. The work of J.C. and B.F.G. was 
supported by the U.S. Department of Energy under contract W-7405-ENG-36, and the work 
of R.S. was supported by DOE contract DE-AC05-84ER40150 under which the Southeastern 
Universities Research Association (SURA) operates the Thomas Jefferson National Accel- 
erator Facility. Finally, some of the calculations were made possible by grants of computing 
time from the National Energy Research Supercomputer Center. 



17 



REFERENCES 



[1] A.R. Berdoz et al, |aucl-ex/0107014| , submitted to Phys. Rev. Lett. (2001). 
[2] W. M. Snow, et al, Nucl. Instrum. Meth. A440, 729 (2000). 

[3] B. Wojtsekhowski and W.T.H. van Oers (spokespersons), JLAB letter-of-intent 00-002. 

[4] R. Hasty et al, Science 290, 2021 (2000). 

[5] D.T. Spayde et al, Phys. Rev. Lett. 84, 1106 (2000). 

[6] D.E. Nagle et al, in High Energy Physics with Polarized Beams and Polarized Targets- 
1978, edited by G.H. Thomas, AIP Conference Proceedings No. 51 (American Institute 
of Physics, New York, 1979), p. 218. 
[7] R. Balzer et al, Phys. Rev. C 30, 1409 (1984). 
[8] P.D. Eversheim et al, Phys. Lett. B 256, 11 (1991). 
[9] S. Kistryn et al, Phys. Rev. Lett. 58, 1616 (1987). 
[10] V. Yuan et al, Phys. Rev. Lett. 57, 1680 (19886). 
[11] M. Simonius, Phys. Lett. 41B, 415 (1972). 

[12] V.R. Brown, E.M. Henley, and F.R. Krejs, Phys. Rev. C 9, 935 (1973); Phys. Rev. Lett. 
30, 770 (1973). 

[13] E.M. Henley and F.R. Krejs, Phys. Rev. D 11, 605 (1975). 

[14] T. Oka, Prog. Theor. Phys. 66, 977 (1981). 

[15] D.E. Driscoll and G.A. Miller, Phys. Rev. C 39, 1951 (1989). 

[16] B. Desplanques, J.F. Donoghue, and B.R. Holstein, Ann. Phys. (NY.) 124, 449 (1980). 
[17] R.B. Wiringa, V.G.J. Stoks, and R. Schiavilla, Phys. Rev. C 51, 38 (1995). 
[18] V.G.J. Stoks, R.A.M. Klomp, C.P.F. Terheggen, and J.J. de Swart, Phys. Rev. C 49, 
2950 (1994). 

[19] R. Machleidt, Phys. Rev. C 63, 024001 (2001). 

[20] J.R. Bergervoet et al, Phys. Rev. C 41, 1435 (1990). 

[21] V.G.J. Stoks, R.A.M. Klomp, M.C.M. Rentmeester, and J.J. de Swart, Phys. Rev. C 

48, 792 (1993). 
[22] B.v. Przewoski et al, Phys. Rev. C 58, 1897 (1998). 
[23] J.L. Friar, Ann. Phys. (N.Y.) 104, 380 (1977). 
[24] J.L. Forest, Phys. Rev. C 61, 034007 (2000). 
[25] S.A. Coon and J.L. Friar, Phys. Rev. C 34, 1060 (1986). 
[26] R.Schiavilla, Nucl. Phys. A 689, 84c (2001). 

[27] M.L. Goldberger and K.M. Watson, Collision Theory (Wiley, New York, 1964). 
[28] M. Abramowitz and I.S. Stegun, Handbook of Mathematical Functions (Dover, New 
York, 1974). 

[29] J.T. Holdeman and R.M. Thaler, Phys. Rev. Lett. 14, 81 (1965). 
[30] CM. Vincent and S.C. Phatak, Phys. Rev. C 10, 391 (1974). 

[31] W. Glockle, The Quantum Mechanical Few-Body Problem ( Springer- Verlag, Berlin, 
1983). 

[32] R. Machleidt, private communication. 
[33] W.T.H. van Oers, private communication. 
[34] M. Simonius, Can. J. Phys. 66, 548 (1988). 
[35] R.V Reid, Ann. Phys. (N.Y.) 50, 411 (1968). 
[36] R.B. Wiringa, private communication. 

[37] R. Machleidt, K. Holinde, and Ch. Elster, Phys. Rep. 149, 1 (1987). 



18 



[38] R.A. Arndt, J.S. Hyslop III, and L.D. Roper, Phys. Re v. D 35, 128 (198 7). 
[39] W.C. Haxton, CP. Liu, M.J. Ramsey-Musolf, preprint |nucl-th/ 0109011 
[40] J. Lang et al, Phys. Rev. C 34, 1545 (1986). 
[41] B. Heckel, private communication. 



19 



TABLES 



TABLE I. Values used for the strong- and weak-interaction coupling constants of the p- and 
cj-meson to the nucleon, see text. 







10 7 /ig p (DDH-adj) 


I0 7 h p a p (DDH-orig) 


A Q (GeV/c) 


p 0.84 
u 20. 


6.1 
0. 




-22.3 
+5.17 


-15.5 
-3.04 


1.31 
1.50 


TABLE II. Labeling of channels. 


J 




1 




a 
2 


3 



1 
2 
3 
4 




% 
3 Pi 
3 P 2 
3 F 3 
3 F 4 




3 Po 
3 F 2 
3 H 4 




TABLE III. 

respectively. 


Classification of channel mixings for J 


even: PC or PV if induced by 


PP PV 

v or v , 


J 




12 




coupling 
13 


23 



2 
4 




PV 
PC 
PC 
PC 




PV 
PV 
PV 


PV 
PV 
PV 
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FIG. 1. Total nuclear asymmetries obtained with various combinations of strong- 
weak-interaction potentials are compared to data, see text. 
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FIG. 2. Contributions to the total nuclear asymmetry obtained by including only the J 
channel, and by adding the J=2 channels, and finally all even J-channels up to J m ax=8. 
AV18- and DDH-adj potential combination is used, thin solid line in Fig. |l|. 
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FIG. 3. Sensitivity of the total nuclear asymmetry to modifications of the cutoff-parameters 
Ap and A w in the DDH potential. Both cut-offs are multiplied by i? C ut 5 see text. For each case the 
couplings are then adjusted to approximately reproduce the results obtained with the AV18 and 
DDH-adj potential combination (corresponding to R cu t=l) in Fig. |T[ 
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FIG. 4. The total nuclear asymmetry, see text, is compared to the total asymmetry obtained 
by ignoring the Coulomb potential. The AV18 and DDH-adj potential combination is used, thin 
solid line in Fig. |I[ 
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FIG. 5. Angular distributions obtained with (curve labeled C) and without (curve labeled no-C) 
inclusion of the Coulomb amplitude at 45 MeV and 221 MeV, see text. The AV18 and DDH-adj 
potential combination is used, thin solid line in Fig. [l]. 
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FIG. 6. The total nuclear asymmetry (thick solid line) is compared with the total asymmetries 
measured in transmission experiments with critical angles 9o=2°, 5°, and 10° (curves labeled A |>2, 
A |>5, and A |>io, respectively). The AV18 and DDH-adj potential combination is used, thin solid 
line in Fig. |l|. 
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FIG. 7. Angular distributions of the longitudinal asymmetry obtained with either the AV18 or 
BONN potential in combination with the DDH-adj potential. 
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FIG. 8. Curves of constant total x 2 obtained by analyzing the experimental pp data with the 
AV18 model, and p- and u>meson strong-interaction couplings in the DDH potential from the 
Bonn-2000 model. The curves indicate surfaces of total \ 2 = 1)2,3,4, and 5 for various values of 
h p p p and hPJ>. 
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